Defect Formation Rates in Cosmological First-Order Phase Transitions 
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In cosmological first-order phase transitions, the progress of true-vacuum bubbles is expected to 
be significantly retarded by the interaction between the bubble wall and the hot plasma. It has 
been claimed that this leads to a significant reduction in the number of topological defects formed 
per bubble, as a result of phase equilibration between bubbles. This claim has been verified for 
spontaneously-broken global symmetries. We perform a series of simulations of complete phase 
transitions in the 2 -I- 1-dimensional i7(l)-Abelian Higgs model, for a range of bubble wall velocities, 
in order to obtain a quantitative measure of the effect of bubble wall speed on the number density 
of topological defects. We find that the number of defects formed is i) significantly lower in the 
local than the global case and ii) decreases exponentially as a function of wall velocity. Slow-moving 
bubbles also lead, however, to the nucleation of more bubbles before the phase transition is complete. 
Our simulations show that this is in fact the dominant effect, and so we predict more defects per 
unit volume as a result of the sub-luminal bubble wall terminal velocity. 
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I. INTRODUCTION 

As the early Universe expanded and cooled, it is 
expected to have undergone a series of phase transi- 
tions at each of which a symmetry which had been 
thermally-restored became spontaneously-broken. At 
each of these phase transitions, depending on the nature 
of the symmetry-breaking involved, topological defects 
may or may not have formed [Q. Phase transitions are 
labelled first- or second-order according to whether the 
field-space location of the true vacuum state changes dis- 
continuously or continuously as the critical temperature 
is crossed. First-order phase transitions, with which this 
work is concerned, proceed by bubble nucleation and ex- 
pansion. Depending on the details of the theory, when 
at least (4 — n) bubbles collide, an n— dimensional topo- 
logical defect may form in the region between them. In 
this paper, we simulate a series of complete cosmologi- 
cal first-order phase transitions in order to quantify the 
initial density of topological defects. 

We take as our model the simplest spontaneously- 
broken gauge symmetry: the Abelian Higgs model, which 
has a local U{\) symmetry, with Lagrangian 
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where 
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The detailed form of the effective potential will 
depend upon the particular particle-physics model being 
considered, but in order to be able to study the generic 
features of a first-order phase transition we shall take V , 
following 101 and IHl, to be 
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where in a realistic model, e — e{T) (x {Tc — T). The 
potential V has a local minimum false-vacuum state 
at $ = which is invariant under the U{\) symme- 
try, and global minima true-vacuum states on the circle 
|$| = ptv = ('7/4) (3 -I- £ -f Vl + 6e -|- e^) which possess 
no symmetry. The dimensionless parameter e is respon- 
sible for lifting the degeneracy between the two sets of 
minima - the greater e, the greater the potential dif- 
ference between the false- and true-vacuum states, and 
hence the faster the bubbles will accelerate. This model 
admits string-solutions, or in the 2 -I- 1-dimensional case 
considered here, vortices. 

By making the field and coordinate transformations 
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it is possible to set A and rj to unity, so that the potential 
is parametrised only by e, and hereafter we shall use these 
transformed variables. 

The bubble nucleation rate per unit time per unit vol- 
ume and the field-space profile of the bubble wall are 
obtained from the 'bounce' (i.e. least action) solution of 
the Euclidean action Q|. In rf+l-dimensions the modulus 
p of the field <i> is given by the solution of the equation 
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with boundary conditions 
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r=0 

Ignoring quantum fluctuations, the phase 9 is constant 
within each bubble, and uncorrelated between spatiaUy- 
separated bubbles. Any non-zero gauge fields in the nu- 
cleation configuration will make a contribution to the ac- 
tion and hence the nucleation of bubbles with non-zero 
gauge fields is exponentially suppressed. When three or 
more bubbles collide, a phase-winding of 2?™ can occur 
around a point, which by continuity must then be at 
<i> = 0. In three spatial dimensions, this topologically- 
stable region of high-energy false vacuum is string-like 
- a cosmic string. If the phases of the bubbles were 
uniformly- distributed, it is simple to calculate that at ev- 
ery three-bubble collision, there is a probability p — 0.25 
that the arrangement of phases would be such as to gen- 
erate the winding required to form a defect. We will 
use this value of p = 0.25 as our benchmark through- 
out this work. This description of defect formation and 
the estimation of the initial defect density is known as 
the Kibble mechanism, and the purpose of this work is 
to investigate to what extent its predictions are modified 
by the early-Universe environment in which cosmological 
phase transitions take place. 

Although two-dimensional (or domain walls) and 
pointlike defects (monopoles) [Hj^ would dominate the 
energy density of the Universe, one-dimensional 'cosmic 
strings' are expected to reach a scaling solution (where 
the ratio of the string energy density to that of the back- 
ground remains constant), and thus need not be patho- 
logical. For many years, cosmic string-seeded perturba- 
tions were considered a viable candidate for the primor- 
dial density fluctuations which led to the formation of 
structure in the Universe. The failure of cosmic string- 
induced CMB anisotropy spectra ||] to match observa- 
tions however, has laid rest to this theory. This 
does not mean though that cosmic strings did not form, 
or that if they did, they are cosmologically-irrelevant. 

Cosmic strings may be either closed loops, or infinitely- 
long strings (see jllj or for cosmic string reviews). 
Closed loops will decay via gravitational radiation, form- 
ing a gravitational-wave background which could be de- 
tectable in the planned experiments LIGO and LISA. 
Cosmic strings have been cited as being responsible for 
baryogenesis [[l3yi4[ , if they trap a Grand Unified gauge 
field with baryon-number violating interactions in their 
core. They could cause the highest energy cosmic rays 
I psf , whose existence has largely defied explanation ~ 
standard methods of ray production such as quasars do 
not appear to be able to generate such events. When 
string loops double back on themselves, forming a 'cusp', 
a jet of ultra-high energy accelerated particles is emitted, 
which could explain the origin of Gamma Ray Bursts p6| . 
It has even been proposed that cosmic strings may consti- 
tute the dark matter - a network of non-Abelian (in 



order to prevent a scaling solution being reached) strings 
would have negative pressure and could therefore explain 
the why the Universe appears to be accelerating. 

In order to be able to assess the significance of cosmic 
strings in the evolution of the early Universe, it is impor- 
tant to be able to estimate the initial defect density accu- 
rately. If the phase transition at which they formed was 
first-order, this would depend on how the phases between 
two or more bubbles interpolate after collision. In partic- 
ular, although strings are in general formed when three or 
more bubbles collide, a simultaneous three-bubble colli- 
sion is unlikely - one would expect in general two-bubble 
collisions, with a third, or fourth bubble colliding some 
finite time later. If the phase inside a two-bubble col- 
lision is able to equilibrate quickly, and before a third 
bubble arrives, there may be a strong suppression of the 
initial string density. The effect of phase equilibration on 
the initial defect density was first investigated by Melfo 
and Perivolaropoulos [Q . They found a decrease of less 
than 10% (i.e. 0.22 < p < 0.25), in models which pos- 
sess a global symmetry and with bubbles which, upon 
nucleation, accelerated up to the speed of light. 

The above description of defect formation, however, 
ignores any effect that the hot-plasma background may 
have on the evolution of the Higgs field, which may be sig- 
nificant in the early Universe. Real-time simulations |l9|] 
and analytic calculations ||2^ for the (Standard Model) 
electroweak phase transition predicted that the bubble 
wall would reach a terminal velocity vtcr 0.1c. Recent 
calculations for the Minimal Supersymmetric Standard 
Model, where there are many more particles, and hence 
more interactions and hence a more viscous plasma, 
show that the bubble wall velocity could be as low as 
Vtcr ~ lO^'^c The reason for this is simple: outside 
the bubble, where the {SU{2) x U{1)) symmetry remains 
unbroken, all fields coupled to the Higgs are massless, ac- 
quiring their mass from the vacuum expectation value of 
the Higgs in the spontaneously-broken symmetry phase 
inside the bubble. Particles outside the bubble without 
enough energy to become massive inside bounce off the 
bubble wall, retarding its progress through the plasma. 
The faster the bubble is moving, the greater the momen- 
tum transfer in each collision, and hence the stronger the 
retarding force. Thus a force proportional to the bubble- 
wall velocity appears in the effective equations of motion. 

Ferrera and Melfo |^ studied bubble collisions in such 
an environment, for theories which possessed a global 
symmetry, and found that decaying phase oscillations 
occur inside a two-bubble collision, which was claimed 
would lead to a suppression of the defect formation rate. 
Kibble and Vilcnkin |22j studied phase dynamics in col- 
lisions of undamped bubbles in models with a local sym- 
metry, and found, analytically, a different kind of decay- 
ing phase oscillation. When the finite conductivity of the 
plasma was included, these oscillations were found not to 
occur. Davis and Lilley extended this work, to include 
slow- moving bubbles in local-symmetry models. Simu- 
lating the collision and merging of two bubbles, it was 
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found that the phase oscillations were suppressed - the 
(gauge-invariant) phase difference between two bubbles 
equilibrated very rapidly upon collision (and even more 
rapidly when the conductivity of the plasma was taken 
into account), with the implication that there would be a 
suppression of defect formation in slow-moving bubbles 
in gauge theories as well. This claim was illustrated by 
an example of an undamped 3-bubble collision which led 
to the formation of a defect, but when identical initial 
conditions were evolved in a damped environment, no 
defect formed. 

Ferrera ||2^] confirmed the hypothesis of ||^ by perform- 
ing a series of simulations of complete global-symmetry 
phase transitions at different values of the bubble wall 
terminal velocity. However, quantitative analysis of the 
defect formation rate in the most realistic scenario cos- 
mologically - a gauge-theory phase transition where the 
bubbles are slowed significantly by the plasma (as might 
be expected at the electroweak- or GUT-scales) - has not 
been studied. 

In this paper we present the results of our investiga- 
tions into the effect of the sub-luminal bubble wall veloc- 
ity in gauge theory phase transitions. We have simulated 
a series of complete phase transitions, from the nucleation 
of the first bubble in an entirely symmetric background 
until > 95% of the available space having been converted 
into the broken-symmetry phase. We have done this (as 
in p3|] , but in a local rather than global symmetry) for 
a range of values of the bubble wall velocity in order to 
obtain a qualitative measure of the dependence of on 
ftci- Our results (which are displayed in Figure]^) show 
that the two hypotheses of ^ do hold: for any value of 
the bubble wall speed there are fewer defects in gauge 
theory phase transitions than in global theory ones, and 
the number of defects formed per bubble decreases with 
decreasing wall velocity. For values of the wall velocity 
which are high compared with those calculated for cosmo- 
logical phase transitions, we find already almost an order 
of magnitude fewer defects. We also confirm the gauge- 
theory analogue of the result of - if the bubbles move 
at the speed of light, there is no significant reduction in 
the defect formation rate due to phase equilibration. 

As well as the possibility of phase equilibration and 
the suppression of defect formation however, slow- moving 
bubbles will also require a larger number of bubbles to 
be nucleated in order to complete the phase transition. 
These two effects act in opposition, and our simulations 
show, for the range of wall velocities considered, that the 
former is dominant, i.e. the net effect of slow-moving 
bubbles is that more defects per unit volume will be 
formed than one would expect were the bubbles to move 
at the speed of light. 

We should note in passing that we have ignored the 
effect of the expansion of the Universe in our work. This 
is a good approximation for phase transitions which take 
place at late times, like the electroweak phase transition. 
At phase transitions which occur earlier, however, the 
Hubble expansion may have a significant effect on bubble 



and phase dynamics. This topic deserves consideration 
on its own. 



II. SIMULATING A PHASE TRANSITION 

In the previous section we have described how the sub- 
luminal terminal velocity of the bubble walls causes the 
gauge-invariant phase difference between two bubbles to 
equilibrate more quickly than in the undamped case. It 
has been demonstrated, qualitatively, how this can lead 
to the suppression of topological defect formation |^ , and 
we have explained how this could be cosmologically sig- 
nificant. In order to set this work on firmer footing, and 
be in a position to assess the cosmological consequences, 
we would like now to be able to quantify these effects. 

In particular, we would like to be able to investigate 
two of the claims made in Q : 

1. For a given wall velocity, fewer topological defects 
are formed in local-symmetry models than in those 
with a global symmetry. 

2. In models where a local symmetry is spontaneously 
broken, the average number of defects formed per 
bubble nucleated is a decreasing function of the 
bubble wall velocity. 

In order then to quantify the effect on defect formation 
probabilities of the slow-moving bubble walls, we must 
perform many simulations of 'complete' phase transitions 
each involving the nucleation of many bubbles. Such an 
investigation was carried out for values of the terminal 
velocity in the range uter = c to Wtcr ~ 0.05c. This paper 
describes the procedure followed and presents the results 
obtained. 

In the standard computer simulations of defect 
formation, causally-disconnected points on the spatial 
lattice are assigned randomly-generated relative phases 
(corresponding to the centres of true-vacuum bubbles in 
a first-order phase transition, or of domains in a second- 
order transition). The phase between sites is then taken 
to vary on the shortest path on the vacuum manifold - 
the so-called geodesic rule. Topological defects will then 
form wherever this geodesic interpolation between sites 
generates a topologically-nontrivial path in the vacuum 
manifold. For a first-order transition this formalism cor- 
responds to true vacuum bubbles nucleating simultane- 
ously, equidistant from all their nearest neighbors. Con- 
sequently all collisions between neighboring bubbles oc- 
cur simultaneously, and the associated phase differences 
are simply given by the differences in the initial assigned 
phases. 

We, on the other hand, wish to investigate precisely the 
effect on the defect formation probability of phase equi- 
libration - something which is only possible as a result 
of the non-simultaneity of three-bubble collisions. 
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The procedure described above was formulated as a 
consequence of the hmitations of computing power avail- 
able at the time. Now, however, with the increased pro- 
cessor capabilities, we have the opportunity to perform 
much more computationally-intensive simulations. As a 
result, instead of following the standard procedure, we 
have used a spatial lattice several times smaller than the 
width of the bubble wall, and evolve the discretised field 
equations in their entirety. This is necessary as it is the 
dynamics of the phase and gauge fields which lead to the 
suppression of defect formation seen in ||] . 



III. PHASE TRANSITION ALGORITHM 



We wish to simulate a first-order phase transition in 
the Abelian-Higgs Model, which possesses a J7(l) gauge 
symmetry and has equations of motion 



p" -{^^e-eA^fp = - 



^V 



^ ^' ^ dp 
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i; - AJ' -d,(d-A)^ -2ep'd,e, 
where we will take the potential V to be given by 
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with A and 77 re-scaled to unity in the manner described in 
the Introduction. As we have stated however, we would 
like to investigate the behaviour of the phase in collisions 
of slow-moving bubbles. For a given theory, by consider- 
ing the Boltzmann equations for scattering off the Higgs 
field, it is possible to calculate the terminal velocity of 
the bubble wall Since we are not concerned here 

with the parameters of a specific particle-physics model, 
we choose instead to use a single damping parameter T 
to model the interaction of the Higgs with the plasma. 
In the introduction we claimed that the plasma would 
introduce a term proportional to the bubble-wall veloc- 
ity into the equations of motion. Since the phase of 
the Higgs field is not affected by the effects described, we 
assume that the plasma couples only to the modulus p. 
We then have an effective equation of motion for p 

BV 

p-p" + rp-{d,0-eA,fp = - — , (13) 

replacing equation (^. Equations (|l^) and ( pj] ) remain 
unchanged. 

A damping term of this form has been used by several 
authors , , , and has also been derived from the 
stress-energy of the Higgs, assuming a coupling to the 
plasma [Q. Heckler |Q estimates F ~ S^f^c for the 
electroweak phase transition, by comparing the energy 
generated by the frictional damping with the pressure on 
the wall due to the damping. 



FIG. 1. Log-log plot of the number of defects formed per 
nucleated bubble, rid, vs the friction coefficient gamma (from 
[psi). The two solid lines correspond to the least squares best 
fits for the two different regimes Ud oc 7~"i-^. The "low" 
friction regime has a slope of qi — 0.53 ±0.05, whereas in the 
high friction regime Ud decays with 02 ~ 0.8 ± 0.3. 

The effect of this damping term is that instead of ac- 
celerating up to the speed of light, the bubble walls reach 
a terminal velocity Vtcr < c. An exact relation between 
F and Wter can be calculated Q|, but the approximation 
Wjgj. ~ 1/F is sufficiently accurate for our purposes. 

The results of Ferrera |^ are shown in Figure |^ - 
the number of defects formed per bubble nucleated, Ud, 
decreases exponentially with decreasing bubble wall ve- 
locity. For terminal velocities vtcr ~ 0.01c, a decrease of 
over an order of magnitude in the defect formation rate 
(i.e. p < 0.03) was measured. 

We will then follow a modified version of Ferrera's pro- 
gram (as in his work using only 2-1-1 dimensions, again 
because of computational constraints - the gauge the- 
ory has many more degrees of freedom than the global 
theory and thus requires significantly more memory and 
computer time) of nucleating bubbles at random points 
in space and time during the course of the simulation. 
Our algorithm is as follows: 

1. for the specific parameters of the potential under 
investigation, compute the bubble wall profile (nu- 
merically) from equation (^. 

2. generate a population of time ordered bubble nu- 
cleation events distributed randomly within some 
finite volume of 2-1-1 dimensional space-time, as- 
signing a random phase to each bubble. 

3. start with an initial state for which the field is in 
the false vacuum (i.c $ = 0) in all the simulation 
volume, and, assuming that there is no primordial 
magnetic field, the gauge fields are everywhere zero. 
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4. at the beginning of each time step check whether 
there are any bubbles to be nucleated at that time: 

(a) if there are any bubble nucleation events, nu- 
cleate only those bubbles that would fall en- 
tirely in a false vacuum region and discard the 
rest (i.e., avoid superimposing new bubbles on 
regions that are already in the true vacuum). 

5. evolve the resulting field configuration to the next 
time step, following the field equations (|l^), ( [lO| ) 
and (|l^) (discretised in the gauge-invariant way de- 
scribed in, for example Q), using a Runge-Kutta 
fourth-order algorithm. 

In order to minimise boundary effects we nucleated 
bubbles only within a box which was a few bubble radii 
smaller than our simulation volume, as was done by Sri- 
vastava , and also used this smaller box to determine 
when the phase transition was complete. 

To generate the bubble nucleation events, we first fix 
their number, then choose at random the space-time 
points at which they take place within the simulation 
volume. Since we only nucleate those bubbles that fall 
within the false vacuum and discard the rest, at later 
times into the transition it will become increasingly dif- 
ficult for the new bubbles to find themselves in the false 
vacuum, and consequently fewer and fewer bubbles will 
be nucleated. Also, and although in a realistic situation 
one would expect the nucleation rate to vary with the 
amount of dissipation present in the system, we kept the 
nucleation rate constant throughout the simulation se- 
ries. This is because we are concerned primarily with 
investigating the number of vortices produced per nu- 
cleated bubble, Ud - since computational resources are 
finite, we have chosen not to explore the effect of varying 
the nucleation rate at this stage. 

The only constraint on the nucleation rate is that we 
wish to ensure that on average the bubbles reach their 
terminal velocity before colliding. The bubbles collide, 
on average, when their radius is of the order of the initial 
spatial separation of nucleation events 



(14) 



where d is the number of spatial dimensions. 

We can parametrise the radius at a given time by three 
quantities: the initial radius tq, the terminal velocity Wtcr 
and the time taken to reach terminal velocity iter- These 
three parameters are easily obtained from simulations. 
Assuming that the bubble wall accelerates uniformly up 
to wtcr and then moves with constant velocity, we get 



r{t) ro + Wtcr t 



^tcr 

Then the condition <coU > iter can be expressed 
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and so the value of used was chosen to satisfy this 
constraint for all values of T. 

Similarly, we have chosen not to investigate the effect 
of crossing the Bogomol'nyi limit - at the critical cou- 
pling /3 — A/2e^ ~ P^^'' of vortices remain stationary 
(see, for instance |1^). Above this limit they repel one 
another, and below they attract. We have chosen, for 
reasons described in section EV, parameters correspond- 
ing to /3 < 1 for our simulations. 



IV. CHOICE OF SIMULATION PARAMETERS 

The intention of this project is to investigate the ef- 
fect of bubble terminal velocity on the defect formation 
process. It is therefore necessary to ensure that the bub- 
bles, once nucleated, reach their terminal velocity t;tor as 
rapidly as possible in order to decrease the probability of 
bubbles colliding at speeds lower than vta and contam- 
inating our results. This is achieved by choosing a rel- 
atively high value of the asymmetry parameter, e — 0.8 
- Figure |^ shows the effect of e on the acceleration of a 
bubble. This choice also has two advantages: the ratio 
of the radius of a nucleated critical bubble to its thick- 
ness is relatively small (compared to lower values of e), 
which aids numerical accuracy (by decreasing the spa- 
tial resolution required to produce stable results) and it 
is the value chosen by Ferrera for his simulations, thus 
enabling a direct comparison of the results in global and 
local-symmetry phase transitions. Figure |^ shows the 
difference in the bubble wall profile for small and large 
values of e. 

For the gauge couphng constant, we took e = 0.25 
throughout. At each value of F a number of phase tran- 
sitions at e = (i.e a global symmetry model) were sim- 
ulated to ensure consistency with the results of Ferrera 
(Figure |^). Since we primarily wish to investigate the 
effect of the bubble wall speed on the number of defects 
formed, we have chosen to perform all our calculations 
with the conductivity cr = 0. 



V. ENDING THE PHASE TRANSITION 

Vortices have dynamics, and as such it is important 
that we specify when we intend to count them. If the 
simulation were allowed to run indefinitely, all vortices 
formed would eventually annihilate, leaving none - mag- 
netic fiux is conserved, and since we started out with zero 
total flux, at any given time there must be an equal num- 
ber of vortices and antivortices. Following j23|, we have 
decided to terminate the simulation when 95% of the sim- 
ulation volume had been converted to true vacuum. As 
a check, the number of vortices and antivortices present 
was calculated at regular intervals and stored to ensure 
that a) the number of defects was an increasing func- 
tion of time, and b) the number of defects had remained 
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Time 

FIG. 2. Illustration of effect of asymmetry parameter e on 
the acceleration of the bubble wall. The solid line shows the 
evolution of wall speed with time for e = 0.4. This bub- 
ble, with a small radius and large difference in energy density 
across its membrane, accelerates rapidly towards the speed of 
light. The dashed line shows the wall speed for e = 0.1. The 
speed of this bubble wall approaches unity much more slowly. 




FIG. 3. Comparison of wall profiles for different parame- 
ters. The dashed line shows the profile for e — 0.1, which has 
a thin wall and large radius. The solid line shows the profile 
for e — 0.4 and the dotted line that for e = 0.8, the value 
used in this paper. The radii of the bubbles decreases, and 
the thickness of the wall increases with increasing e, making 
high values of e most suitable for these simulations. 



constant for some time before the phase transition was 
stopped. Any uncertainty introduced by this choice of 
95% completion would have been compensated by the 
large number of simulations performed for each value of 
r. Indeed a glance at the cumulative vortex count shows 
that the number of defects present in the simulation vol- 
ume does not change significantly if the truncation time 
is increased or decreased by 10%. 

At the end of the phase transition, the number of vor- 
tices and antivortices produced were counted, using an al- 
gorithm which went across the grid checking at each point 
for gauge-invariant phase windings of 27rn around the 
smallest square going clockwise from that point with all 
four corners in true vacuum. That the number of vortices 
and antivortices must be equal is another check on the 
accuracy of the code. 



VI. DETAILS OF COMPUTATION 

The first-order phase transition was simulated for 7 
values of F. For F = 0,2,3,5 and 7.5, 30 simulations 
were performed at each value, whereas for F = 10, only 
5 were performed, due to the increasing amount of time 
taken to complete the phase transition. At each value 
of F, the simulations were run on a smaller and smaller 
lattice spacing with identical initial conditions, until no 
divergence in the results could be detected. 

A lattice spacing of 0.02 was determined to be suffi- 
cient for all values of F and so was used throughout. For 
reasons of simulation time, rather than memory, the spa- 
tial lattice size was then set at 3000^ grid points, giving 
a total physical area of 60^. This box size corresponds to 
around 20 bubble radii, for the case where e — 0.8. The 
time step At was set so that it satisfied the Courant con- 
dition (see, for instance p. 829 of for further details) 
for numerical stability - namely 



At ^ 1 



(17) 



where d is the number of spatial dimensions. Altogether, 
we estimate that this project required around 1300 hours 
of combined CPU time on the supercomputer. 



VII. SIMULATION RESULTS 

The main result of our simulations can be seen in Fig- 
ure ||. Our conclusions are twofold and can be stated 
thus - the number of topological defects formed per bub- 
ble nucleated did decrease significantly as the terminal 
velocity of the wall decreased; the number of topologi- 
cal defects formed per bubble nucleated was significantly 
lower for gauge theories than for global theories. Both of 
these statements support the hypotheses made in 13] . 
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FIG. 4. Snapshots of a phase transition with F = at different times. From top left, they are ai t — 2.2, 3.67, 5.13, 7.33, and 
the phase transition was declared complete by the program ai t — 7.80. The length of the arrows is proportional to the value 
of |$| = p, and their direction gives the phase 9 of the Higgs field. In particular, the false vacuum (which is $ = 0) corresponds 
to an absence of lines. The first image is very early in the phase transition, when only very few bubbles have been nucleated. 
These bubbles expand, alongside new bubbles which continue to be nucleated. In the second and third of these images, the 
existence at a given instant of bubbles of differing sizes is evident. The final image shows the phase transition almost complete. 



We have shown that the number of defects formed per 
bubble nucleated, rid-, decreases exponentially with bub- 
ble wall velocity (since wtcr ^ 1/r). The cosmologically- 
relevant quantity would be, of course, the number of de- 
fects per unit volume, riy. These are related by 



^ number of bubbles per unit volume r ' 

where r is the average bubble radius on collision. Since 
the only parameters in this system are the damping co- 
efficient F and the bubble nucleation rate per unit time 
per unit volume N, we have 

r^lNT]-^, (19) 

where d is the number of spatial dimensions. If, as we 
have discovered, ~ F^", we have 

n„ - F(3fT-"). (20) 

If we assume a fixed nucleation rate, slow-moving bub- 
bles have two effects with respect to the formation of 
topological defects. The first is that more bubbles must 
be nucleated to complete the phase transition, and the 
second is that phase equilibration suppresses the forma- 
tion of defects in some cases. Our results indicate, since 
we have found a < 1 for d — 2, that the increase in 
the number of bubbles nucleated outstrips the decrease 
in the defect-formation rate due to phase equilibration, 
and that slow-moving bubbles lead to more defects in any 
given volume. 

Here we must point out that this statement is only 
valid for the range of wall velocities we have considered 
(0.05 < Vtor/c < 1). In this range, we have measured 
the exponent aiocai — 0.65 ± 0.04, whereas in the global 
case [HI agiobai = 0.53 ± 0.05 was found for the corre- 
sponding region. For lower terminal velocities, however, 
another regime was discovered, with Ogiobai = 0.8 ± 0.3. 
The crossover between the two regimes corresponded to 
a bubble velocity low enough that the time taken for the 
bubble walls to merge on collision was longer than the 




r 



FIG. 5. Log-log plot of number of defects formed per bub- 
ble nucleated vs. the friction coefficient F. The lower set 
of points (data points marked with a cross) are the results 
of this investigation, the upper set (marked with a circle) are 
those obtained by Ferrera |2^. The two hypotheses of ^ 
are demonstrably supported - namely i) that Ud should de- 
crease with increasing F, and ii) that for a given value of F, 
Ud should be lower for e = 0.25 (the gauge symmetry) than 
for e = the global symmetry. The solid line corresponds to 
the least-squares, best-fit for Ud tx F", with a — —0.65 ±0.04. 

time light took to cross the bubble. We have not been 
able to probe this velocity region, and cannot rule out the 
possibility that for very slowly-moving bubble walls, de- 
fect suppression due to phase equilibration will outstrip 
the effect of more bubbles being nucleated, and so lead 
to fewer defects per unit volume. 

A. Results for F = 

Not shown on Figure ^ (which is drawn with log-log 
axes) is the result for F = ~ for this we obtained 
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Ud = 0.240 ± 0.02. This agrees well with the predic- 
tion of the Kibble mechanism, that on purely geomet- 
rical grounds the number of defects formed per bubble 
should be 0.25. As stated in the Introduction, Melfo and 
Perivolaropoulos found that phase equilibration due 
to non-simultaneous three-bubble collisions did not sig- 
nificantly alter the defect formation rate for global sym- 
metries (e = 0), with bubbles moving at the speed of 

light (r = 0). 

We can now confirm that in our gauge-theory model, 
with bubbles moving at the speed of light, we find no sig- 
nificant deviation from the Kibble mechanism prediction. 
This result, we feel, is significant enough to be stated in 
its own section. 
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VIII. SUMMARY 



In this paper we have given an account of our attempts 
to simulate a first-order phase transition incorporating 
one of the effects of an early-universe environment - a 
retarded bubble wall. We have found, as we predicted, 
that the number of topological defects formed per bubble 
nucleated, n^, decreases with the terminal velocity of the 
bubble wall. Moreover, the number of topological defects 
produced is significantly lower, for every value of F, in the 
gauge theory than in the global theory. 

Slow-moving bubbles allow, in general, more time be- 
tween collisions and so facilitate the suppression of de- 
fect formation due to phase equilibration. Another effect 
of slow-moving bubbles, however, is the need for more 
bubbles to be nucleated in order for the phase transi- 
tion to complete. Our studies show that the surfeit of 
bubbles dominates the suppression of defects, and that 
slow-moving bubbles (0.05 ^ "ytor/c < 1) actually lead to 
more defects in a given volume. 

This investigation is far from complete. We would have 
liked to have been able to probe the region vter < 0.01, as 
this may be more realistic for cosmological scenarios. We 
would also like to have investigated the effect of plasma 
conductivity on the defect formation rates, both sepa- 
rately from and in tandem with the probe on the effect 
of wall terminal velocity. Unfortunately, time and com- 
putational constraints did not allow this up until now. 

Of course, any such cosmological phase transition 
would take place not in a static background, but in an 
expanding universe. The expansion of the Universe has 
not been taken into account here, but could be extremely 
important for any high-energy (e.g. GUT-scale) phase 
transition. This, along with the low-velocity and high- 
conductivity regimes, is something we hope to probe in 
the future. 
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